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Abstract 

This paper introduces the seq_bin meta-constraint 
with a polytime algorithm achieving general- 
ized arc-consistency according to some properties. 
seq_bin can be used for encoding counting con- 
straints such as Change, Smooth or Increas- 
ing_Nvalue. For some of these constraints and 
some of their variants GAC can be enforced with 
a time and space complexity linear in the sum of 
domain sizes, which improves or equals the best 
known results of the literature. 

1 Introduction 

Many constraints are such that a counting variable is equal 
to the number of times a given property is satisfied in a se- 
quence of variables. To represent some of these constraints in 
a generic way, we introduce the SEQ_BiN(iV, X, C, B) meta- 
constraint, where N is an integer variable, X is a sequence of 
integer variables and C and B are two binary constraints. 

Based on the notion C-stretch, a generalization of 
stretch [Pesant, 2001] where the equality constraint is made 
explicit and is replaced by C, seq_bin holds if and only if two 
conditions are both satisfied: (1) TV is equal to the number of 
C-stretches in the sequence X, and (2) B holds on any pair 
of consecutive variables in X. 

Among the constraints that can be expressed thanks to 
seq_bin, many were introduced for solving real-world prob- 
lems, e.g., Change [Cosytec, 1997] (time tabling problems), 
Smooth [Beldiceanu et al, 2010a] (time tabling and schedul- 
ing), or Increasing_N value [Beldiceanu et al, 2010b] 
(symmetry breaking for resource allocation problems). 

The main contribution of this paper is a generic polytime 
filtering algorithm for seq_bin, which achieves generalized 
arc-consistency (GAC) according to some conditions on B 
and C. This algorithm can be seen as a generalization of 
the IncreasingJM value filtering algorithm [Beldiceanu et 
al, 2010b]. Given n the size of X, d the maximum do- 
main size, and Soi the sum of domain sizes, we character- 
ize properties on C and B which lead to a time and space 
complexity in O(Soi). These properties are satisfied when 
seq_bin represents IncreasingJMvalue, and several vari- 
ants of Change (provided its parameter is a monotonic binary 



constraint, e.g., '<', '<', '>', '>'). For these constraints, our 
technique improves or equals the best known results. 

Section 2 provides the definitions used in this paper. Sec- 
tion 3 defines seq_bin and shows how to express well-known 
constraints with seq_bin. Section 4 provides a necessary and 
sufficient condition for achieving GAC. Section 5 details the 
corresponding GAC filtering algorithm. Finally, Section 6 
discusses about related works and Section 7 concludes. 

2 Background 

A Constraint Network is defined by a sequence of variables 
X = [xq,Xi,. . . , x n _i], a sequence of domains V, where 
each D(xi) G T> is the finite set of values that variable Xi 
can take, and a set of constraints C that specifies the allowed 
combinations of values for given subsets of variables, min(x) 
(resp. max(i)) is the minimum (resp. maximum) value of 

D(x). A sequence of variables X 1 = [ Xi , ^2 + 1 , • ■ • ) Xj J, 
< i < j < n — 1 (resp. i > Oori < ?i — 1), is a subsequence 
(resp. a strict subsequence) of X and is denoted by X' C X 
(resp. X 1 C X). A[X] denotes an assignment of values to 
variables in X. Given x G X, A[x] is the value of x in 
A[X] is valid if and only if Vx^ G X, A[xi] G D(xi). An 
instantiation I[X] is a valid assignment of X. Given x G X, 
I[x] is the value of x in I[X}. Given the sequence X and i, 
j two integers such that 0<i<j<n — 1, I[xi, . . . , Xj) 
is the projection of I[X] on [xi, Xj+i, , . . . , Xj}. A constraint 
C(X) G C specifies the allowed combinations of values for 
X. We also use the simple notation C. C(X) defines a subset 
lZc(T>) of the cartesian product of the domains Xl Xi£ xD{xi). 
If X is a pair of variables, then C(X) is binary. We denote by 
vCw a pair of values (v, w) that satisfies a binary constraint 
C. ->C is the opposite of C, that is, ->C defines the relation 
IZ^ci'D) = Ilxiex D(xi) \ IZciTJ). Afeasible instantiation 
I[X] of C\X] is an instantiation which is in lZc(T>). We say 
that I[X] satisfies C{X), or that I[X] is a support on C(X). 
Otherwise, I[X] violates C(X). If C is a binary constraint 
on X = {xi, Xi+i} and v G D(xj) then the set of supports 
such that Xi = v can be considered as a set of values (a subset 
of D(xi + i )). A solution of a constraint network is an instan- 
tiation of all the variables satisfying all the constraints. 

Value v G D(xi), Xi G X, is (generalized) arc-consistent 
(GAC) with respect to C(X) if and only if v belongs to a 
support of C(X). A domain D(xj), Xi G X, is GAC with 
respect to C(X) if and only if Vu G D(xi), v is GAC with 



respect to C(X). C(X) is GAC if and only if Mx t G X, 
D(xi) is GAC with respect to C(X). A constraint network 
is GAC if and only if it is closed for GAC [Bessiere, 2006]: 
Vxi G X all values in -D(xj) that are not GAC with respect 
to a constraint in C have been removed. 

3 The seq_bin Meta- Constraint 

We first generalize the notion of stretches [Pesant, 2001] to 
characterize a sequence of consecutive variables where the 
same binary constraint is satisfied. 

Definition 1 (C-stretch). Let I[X] be an instantiation of the 
variable sequence X = [xo,Xi, . . . , x n _i] and C a binary 
constraint. The C-sequence constraint C(I[X],C) holds if 
and only if: 

• Either n = 1, 

• or n > 1 and Vfc G [0, n — 2] C(I[xk], I[xk+i\) holds. 

A C-stretch of I[X] is a subsequence X' C X such that the 
two following conditions are both satisfied: 

1. The C-sequence C{I[X'\, C) holds, 

2. VAT" such that X' C X" C X the C-sequence 
C{I[X"],C) does not hold. 

The intuition behind Definition 1 is to consider the max- 
imum length subsequences where the binary constraint C is 
satisfied between consecutive variables. Thanks to this gener- 
alized definition of stretches we can now introduce seq_bin. 

Definition 2. The meta-constraint seq_bin(A, X, C, B) is 
defined by a variable N, a sequence of n variables 
X = [xo, X\, . . . , x n -i] and two binary constraints C 
and B. Given an instantiation I\N, Xo, Xi, . . . , x n _i], 
seq_bin(A, X, C, B) is satisfied if and only if for any i E 
[0,n — 2], I[xi] B I[xi+i] holds, and I[N] is equal to the 
number of C-stretches in I[X]. 

The constraint Change was introduced in the context of 
timetabling problems [Cosytec, 1997], in order to put an up- 
per limit on the number of changes of job types during a 
given period. The relation between classical stretches and 
Change was initially stressed in [Hellsten, 2004, page 64]. 
Change is defined on a variable N, a sequence of variables 
X = [xo, xi, . . . , x„_i], and a binary constraint C G {=, / 
,<,>,<,>}. It is satisfied if and only if N is equal to the 
number of times the constraint C holds on consecutive vari- 
ables of X. Without hindering propagation (the constraint 
network is Berge-acyclic), Change can be reformulated as 
seq_bin(AT', X, -iC,true) A [N = N' - 1], where true is 
the universal constraint. 

Smooth(A, X) is a variant of Change(AT, X, C), where 
Xi C Xi+i is defined by \xi — Xj+i| > est, est G N. It is use- 
ful to limit the number of drastic variations on a cumulative 
profile [Beldiceanu et al, 2010a; De Clercq, 2010]. 

As a last example, consider the IncreasingJSIvalue con- 
straint, which is a specialized version of NValue [Pa- 
chet and Roy, 1999]. It was introduced for breaking vari- 
able symmetry in the context of resource allocation prob- 
lems [Beldiceanu et al., 2010b]. Increasing _N value is 
defined on a variable N and on a sequence of variables 



X = [xq, xi, . . . , x n -i). Given an instantiation, Increas- 
ing_N value (TV, X) is satisfied if and only if N is equal to 
the number of distinct values assigned to variables in X, and 
for any i G [0, n — 2], x; < Xj+i. We reformulate Increas- 
ing_N value (Af, X) as seq_bin(A, X, =, <). 

4 Consistency of seq_bin 

We first present how to compute, for any value in a given 
domain of a variable Xi G X, the minimum and maximum 
number of C-stretches within the suffix of X starting at Xi 
(resp. the prefix of X ending at x,) satisfying a chain of bi- 
nary constraints of type B. Then, we introduce several prop- 
erties useful to obtain a feasibility condition for seq_bin, and 
a necessary and sufficient condition for filtering which leads 
to the GAC filtering algorithm presented in Section 5. 

4.1 Computing of the Number of C-stretches 

According to Definition 2, we have to ensure that the chain 
of B constraints are satisfied along the sequence of vari- 
ables X = [xo, xi, . . . , x„_i]. An instantiation I[X] is said 
B-coherent if and only if either n = lorforanyi G [0,n— 2], 
we have I[xi) B I[xi+i]. A value v G D{xi) is said to be 
B-coherent with respect to Xj if and only if it can be part of 
at least one B-coherent instantiation. Then, given an integer 
i G [0, n — 2], if v G D(xi) is .B-coherent with respect to x, 
then there exists w G D(xi+i ) such that v Bw. 

Consequently, within a given domain D(xj), values that 
are not incoherent can be removed since they cannot be part 
of any solution of seq_bin. Our aim is now to compute for 
each £?-coherent value v in the domain of any variable Xj the 
minimum and maximum number of C-stretches on X. 

Notation 1. s(xt,v) (resp. s(xi,v)) is the minimum 
(resp. maximum) number of C-stretches within the se- 
quence of variables [xj, Xj+i, . . . , x n _i] under the hypoth- 
esis that Xi = v. p(xi,v) (resp. p(xi,v)) is the minimum 
(resp. maximum) number of C-stretches within the sequence 
[xo , X\ , . . . , Xi] under the hypothesis that Xi = v. Given 
X = [xo, xi, ... , x n -i], s(X) (resp. s(X)) denotes the min- 
imum (resp. maximum) value of s(xq, v) (resp. s(xq, v)). 

Lemma 1. Given seq_bin(A^ X, C, B) with X = 
[xo,Xi, . . . ,x„_i], assume the domains in X contain only 
B-coherent values. Given i G [0, n — 1] and v G D(xi), 

• If i = n — 1: s(x n -i,v) = 1. 

• Else: 

s(x i ,v)= min ( ^ Bw]A[vCw Mx i+1 ,w)) I \ 

Proof. By induction. From Definition 1, for any v G 
D(x n -i), we have s(x n -i,v) = 1 (i.e., a C-stretch of 
length 1). Consider now Xj G X with i < n — 1, and 
a value v G D(xj). Consider the set of instantiations 
I[xi+i, Xi + 2, . . . , x n -i] that are B-coherent, and that min- 
imize the number of C-stretches in [xi+i,Xi+2, . . . ,x n -i}- 
We denote this minimum number of C-stretches by mins. At 
least one B-coherent instantiation exists since all values in 
the domains of [x,+i, x i+2 , . . . , x n -i] are B-coherent. For 



each such instantiation, let us denote by w the value associ- 
ated with Either there exists such an instantiation 
with mins C-stretches with the conjunction B A C satis- 
fied by (I[xi\, I[xi + i]). Then, s(xi, v) = s(xi + i,w) since 
the first C-stretch of I[xi+i, Xi+ 2 , ■ ■ ■ , %n-x\ is extended 
when augmenting I[xi+i, %i+2, ■ ■ ■ , x n -i] with value v for 
Xi. Or all instantiations I[xi+i, Xi+ 2 , ■ ■ ■ , Xn-i] with mins 
C-stretches are such that C is violated by (/[a;,], I[xi + i]): 
(I[xi\, I[xi+i]) satisfies BA^C. By construction, any instan- 
tiation I[xi, Xi+i, . . . ,x n -i] with I[xi] = v has a number of 
C-stretches strictly greater than mins. Consequently, given 
I[xi+i, Xi+2, ■ ■ ■ , Xn-i] with mins C-stretches, the number 
of C-stretches obtained by augmenting this instantiation with 
value v for Xi is exactly mins + 1. □ 

Lemma 2. Given seq_bin(./V, X, C, B) with X = 
[xo,Xi,... , x n -i], assume the domains in X contain only 
B-coherent values. Given i G [0, n — 1] and v G D(xi): 

• Ifi = n-1: s(x n -i,v) = 1. 

• Else: 

s( Xi ,v)= max ( ^IvBvMvCnMxi+l,™)) \ 

Given a sequence of variables [xq, xi, . . . , x n -i] such that 
their domains contain only B-coherent values, for any Xi in 
the sequence and any v G D{xi), computing p(xi, v) (resp. 
p(xi,v)) is symmetrical to s(xi,v) (resp. s(xi,v)). We sub- 
stitute min by max (resp. max by min), Xi+i by x±-\, and 
vRw by wRv for any R G {B, C, ^C}. 

4.2 Properties on the Number of C-stretches 

This section provides the properties linking the values in 
a domain D(x{) with the minimum and maximum number 
of C-stretches in X. We consider only incoherent values, 
which may be part of a feasible instantiation of seq_bin. Next 
property is a direct consequence of Lemmas 1 and 2. 

Property 1. For any B-coherent value v in D(xi), with re- 
spect to x^ s(xi , v) < s(xi ,v). 

Property 2. Consider seq_bin(./V, X, C, B), a variable x t G 
X (0 < i < n — 1), and two B-coherent values V\,V2 € 
D(xi). If i = n — 1 or if there exists a B-coherent w € 
D(xi+i) such that v\Bw and v^Bw, then s(xi,Vi) + 1 > 
s(x l ,v 2 ). 

Proof. Obviously, if i = n — 1. If v\ = v<x, by Property 1 
the property holds. Otherwise, assume there exist two values 
v\ and t>2 such that 3w G D(xi+\) for which v\Bw and 
viBw, and s(xi,vi) + 1 < s(xi,V2) (hypothesis H). By 
Lemma 2, s(xi , v\ ) > s(xi+i, w). By Lemma 1, s(xj, v 2 ) < 
s(xi_)_i , w) + 1. From hypothesis H, this entails ~s(xi+i , w) + 
1 < s(xi+i , w) + 1, which leads to s(xi+i, w) < s(xi+i , w), 
which is, by Property 1, not possible. □ 

Property 3. Consider seq_bin(./V, X, C, B), a variable x t G 
X (0 < i < n — 1), and two B-coherent values v\,V2 G 
D(xi). If either i = n — 1 or there exists B-coherent 
w G D{xi + \) such that v\Bw and v 2 B w then, for any 
k G [mm(s(xi,Vi),s(xi,v 2 )),max(s(xi,Vi),s(xi,V2))], ei- 
ther k G \s(xi,vi),s(xi,Vi)] orke [s(xi,v 2 ),s(xi,v 2 )]. 



Proof. Obviously, if i = n — 1 or v\ = v 2 . If 
\a(xi,vi),s(xi,Vi)] n \s(xi,v 2 ),s(xi,V2)] is not empty, 
then the property holds. Assume \s(xi, V\), s(xi, 
and [s(xi,v 2 ),s(xi,v 2 )] are disjoint. W.l.o.g., assume 
s(xi,v\) < s(x i} v 2 ). By Property 2, s(xi,Vi) + 1 > 
s(xi,v 2 ), thus s(xi,vi) = s(xi,v 2 ) — 1. Either k G 
\s(xi,vx),s(xi,vi)] or k G [s(x i ,v 2 ),s(x l ,v 2 )] (there is no 
hole in the range formed by the union of these intervals). □ 

4.3 Properties on Binary Constraints 

Property 3 is central for providing a GAC filtering algorithm 
based on the count, for each incoherent value in a domain, of 
the minimum and maximum number of C-stretches in com- 
plete instantiations. Given seq_bin(A^, X, C, B), we focus on 
binary constraints B which guarantee that Property 3 holds. 

Definition 3. [Van Hentenryck et ai, 1992] A binary con- 
straint F is monotonic if and only if there exists a total order- 
ing -< of values in domains such that: for any value v and any 
value w, vFw holds implies v 1 Fw' holds for all valid tuple 
such that v'-<v and w-< w' . 

Binary constraints <, >, < and > are monotonic, as well 
as the universal constraint true. 

Property 4. Consider SEQ_BiN(iV, X, C, B) such that all non 
B-coherent values have been removed from domains of vari- 
ables in X. B is monotonic if and only if for any variable 
xi G X, < i < n — I, for any values Vi,v 2 G D(xi), there 
exists w G D(xi+\) such thatv\Bw and v 2 Bw. 

Proof. (=>) From Definition 3 and since we consider only in- 
coherent values, each value has at least one support on B. 
Moreover, from Definition 3, {w | v 2 Cw} C {w | viCw} or 
{w | V\Cw} C {w | v 2 Cw}. The property holds. (-4=) Sup- 
pose that the second proposition is true (hypothesis H) and 
B is not monotonic. From Definition 3, if B is not mono- 
tonic then 3vi and v 2 in the domain of a variable x- t G A' 
such that, by considering the constraint B on the pair of vari- 
ables (xi,Xi + i), neither {w | v 2 Cw} C {w | ViCw} nor 
{w | ViCw} C {w | v 2 Cw}. Thus, there exists a support 
v\Bw such that (v 2 ,w) is not a support on B, and a support 
v 2 Bw' such that (vi, w') is not a support on B. We can have 
D(xi+i) = {w, w'}, which leads to a contradiction with H. 
The property holds. □ 

4.4 Feasibility 

From Property 4, this section provides an equivalence rela- 
tion between the existence of a solution for seq_bin and the 
current variable domains of X and N. Without loss of gen- 
erality, in this section we consider that all non S-coherent 
values have been removed from domains of variables in X. 
First, Definition 2 entails the following necessary condition 
for feasibility. 

Proposition 1. Given seq_bin(A^, X, C, B), if s(X) > 
max(D(N)) ors(X) < min(D(N)) then seq_bin fails. 

D(N) can be restricted to \s(X),s(X)}, but D(N) may 
have holes or may be strictly included in [s(X),s(X)]. We 
have the following proposition. 



Proposition 2. Consider seq_bin(A, X, C, B) such that B 
is monotonic, with X = [xq, x±, . . . , x n -±]. For any in- 
teger k in \s(X),s(X)] there exists v in D(xq) such that 
k e \s(x ,v),s(x ,v)}. 

Proof. Let V\ <E D(xo) a value such that s(xq, i>i) = s(X). 
Let V2 € D(xq) a value such that s(xq,V2) — s(X). By 
Property 4, either n = 1 or 3w S -D(xi) such that viBw 
and V2BW. Thus, from Property 3, Vfc e [s(X) ,s(X)], either 
fc 6 [s(a;o,«i),s(xo,«i)] or k 6 [s(x , u 2 ), s(a; , u 2 )]- □ 

By Proposition 2, any value for N in D(iV) n \s(X),s(X)} 
is generalized arc-consistent provided a property is satisfied 
on the instance of seq_bin we consider: given a variable x%, 
for any value v in D(xi) and for all k € [s(xi,v),s(xi,v)], 
there exists a solution of seq_bin(A, X, C, B) with exactly 
k C-stretches. 

Definition 4. The constraint seq_bin(7V, X, C, B) is 
counting-continuous if and only if for any instantiation I[X] 
with k C-stretches, for any variable Xi £ X, changing the 
value of Xi in I[X] leads to a number of C-stretches equal 
either to k, or to k + 1, or to k — 1. 

Property 5. Consider seq_bin(A, X, C, B) such that B is 
monotonic, with X = [xq, X\, . . . , x n -i], Xi a variable and 
v G D{xi). If seq_bin(A, X, C, B) is counting-continuous 
then there exists for any integer k € [s(xj, v), s(Xi, v)] an 
instantiation I[{xi, . . . , with exactly k C-stretches. 

Proof. By recurrence, we assume that the property is true 



for all instantiations of [xj , . . 



_i] such that j > 



(the property is obviously true if j = n — 1). At step 

1, we assume that there exists k € [s(xi, v), ~s(xi, v)] such 
that there is no instantiation I[{xi, . . . , x n -i}] with k C- 
stretches, while seq_bin(A, X, C, B) is counting-continuous 
(hypothesis). We prove that this assumption leads to a con- 
tradiction. By Lemmas 1 and 2, there exists an instantia- 
tion I'[{xi, . . . , x n -i}} with s(xi, v) C-streches and I'[xi] = 
v, and there exists an instantiation I"[{xi, . . . , i n _i}] with 
~s(xi,v) C-stretches and /"[x^ = v. Thus, by hypothesis, 
k > s(xi,v) and k < s(xi, v). We have s(xi, v) > s(xi,v) + 

2. By Property 3 and since the property is assumed true for 
all instantiations of [x i+ i , . . . , x n -i), there exists at least one 
pair of values w\, W2 in D(xi+\) such that s(xi+i, W2) = 
s(x i+ i,wi) + 1, (v, wi) satisfies C and (v,W2) violates C 
(this is the only possible configuration leading to the hy- 
pothesis). In this case SEQ_BiN(iV, X, C, B) is not counting- 
continuous: Given an instantiation I[xi, . . . , x n -i] with 
I[xi] = v and /[x^i] = Wi, changing w\ by W2 for Xi + i 
increases the number of C-stretches by 2. □ 

Proposition 3. Given an instance of seq_bin(7V, X, C, B) 
which is counting-continuous and such that B is mono- 
tonic, SEQ_BiN(iV, X, C, B) has a solution if and only if 
\s(X),s(X)}nD(N)^$. 

Proof. (=>) Assume SEQ_BiN(iV, X, C, B) has a solution. Let 
/[{A^} U X] be such a solution. By Lemmas 1 and 2, the 
number of C-stretches I[N] belongs to [s(X),s(X)}. (^=) 
Let k e \s(X),s(X)] n D(N) (not empty). From Proposi- 
tion 2, for any value k in \s(X),s(X)], 3v e D(xq) such 



that k S [s(x ,v),s(x ,v)]. Since seq_bin(A, X, C, B) 
is counting-continuous, there exists an instantiation of X 
with k C-stretches. By Definition 2 and since Lemmas 1 
and 2 consider only incoherent values, there is a solution 
of seq_bin(A, X, C, B) with k C-stretches. □ 

4.5 Necessary and Sufficient Filtering Condition 

Given seq_bin(A, X, C, B), Proposition 3 can be used to fil- 
ter the variable N from variables in X. Propositions 1 and 2 
ensure that every remaining value in [s{X),s(X)] f] D(N) 
is involved in at least one solution satisfying seq_bin. We 
consider now the filtering of variables in X. 

Proposition 4. Given an instance of seq_bin(A, X 1 C, B) 
which is counting-continuous and such that B is monotonic, 
let v be a value in D{xi), i £ [0, n — 1]. The two following 
propositions are equivalent: 

1. v is B-coherent and v is GAC with respect to seq_bin 

\p(xi,v) + s(xi,v) - 1, 



p(xi,v) + s(xi,v) - 1 



n D(N) ^ 



Proof. If v is not incoherent then, by Definition 2, v is not 
GAC. Otherwise, p(xi,v) (resp. s(xi,v)) is the exact min- 
imum number of TJ-stretches among .B-coherent instantia- 
tions I[xq, xi, . . . , Xi] (resp. I[xi, Xi + \, . . . , x„_i]) such that 
I[xi] = v. Thus, by Lemma 1 (and its symmetrical for pre- 
fixes), the exact minimum number of C-stretches among in- 
coherent instantiations I[xq, X\, . . . , x n -i] such that I[xi] = 
v is p(xi,v) + s(xi,v) — 1. Let T)^ ^ C V such that all 
domains in T>u v \ are equal to domains in T> except D{xi) 
which is reduced to {v}. We call Xu v -\ the sequence of 
variables associated with domains in T>u v \. By construction 
p(xi,v) + s(xi 7 v) — 1 = s(X(j^)). By a symmetrical rea- 
soning, p(xi, v) + s(xi, v) — 1 = s(Xf i v \). By Proposition 3, 
the proposition holds. □ 

The "— 1" in expressions p(xi,v) + s(xi,v) — 1 and 
p(xi,v) + s(xi,v) — 1 prevents us from counting twice 
a C-stretch at an extremity x,i of the two sequences 
[x ,xi,...,Xi] and [x i} x i+ \, . . . ,x n -i]. 

5 GAC Filtering Algorithm 

Based on the necessary and sufficient filtering condition of 
Proposition 4, this section provides an implementation of the 
GAC filtering algorithm for a counting-continuous instance 
of seq_bin(A, X, C, B) with a monotonic constraint B. 

If B g {<,>,<,>, true} then the total ordering -< 
entailing monotonicity of B is not the natural order of in- 
tegers. In this case, if -< is not known, it is necessary 
to compute such an ordering with respect to all values 
in Uigfcn-i] (D(xi)), once before the first propagation of 
seq_bin. Consider that the two variables of B can take any 
value in U$ g ro n _i] (D(xi)): Due to the inclusion of sets of 
supports of values (see Definition 3), the order remains the 
same when the domains of the variables constrained by B do 
not contain all values in U ie [ 0i „_i] (D(xi)). 

To compute -<, the following procedure can be used: 
Count the number of supports of each value, in 0(d 2 ) time 



(recall d is the maximum domain size of a variable in X), 
and sort values according to the number of supports, in 

0(| Uiep.n-1] {D{xi))\log{\ U iG[0i „_ 1] (D( Xl ))\)) time. 

Then, given the sequence of variables X, the algorithm is 
decomposed into four phases: 

© Remove all non incoherent values in the domains of X. 

© For all values in the domains of X, compute the mini- 
mum and maximum number of C-stretches of prefixes 
and suffixes. 

® Adjust the minimum and maximum value of N with 
respect to the minimum and maximum number of 
C-stretches of X. 

© Using the result phase © and Proposition 4, prune the 
remaining £>-coherent values. 

With respect to phase ©, recall that B is monotonic: Ac- 
cording to -<, for any pair of variables (xi,Xi+i), 3vq in 
D(xi) such that Vvj G D(x[), vj ^ vq, Vj has a set of 
supports on included in the supports of Vo on 

B(xi,Xi+i). By removing from D(xi+\) non supports of vq 
on B(xi,Xi+i) in 0(\D(xi+i)\), all non incoherent values 
of D(xi+i) with respect to B(xi,Xi+i) are removed. By re- 
peating such a process in the two directions (starting from 
the pair (x n -2, x n -x) and from the pair (xq, all non B- 
-coherent values can be removed from domains in 0(EdO 
time complexity. 

To achieve phase © we use Lemmas 1 and 2 and their sym- 
metrical formulations for prefixes. Without loss of generality, 
we focus on the minimum number of C-stretches of a value Vj 
in the domain of a variable Xi, i < n— 1, thanks to Lemma 1. 
Assume that for all w G D(xi + i), s(xi+i,w) has been com- 
puted. If there is no particular property on C, the supports 
Sj G D(x i+ i) of Vj on C(xi,x i+ i) A B{x i) x i+ i) and the 
subset -iSj G D{xi+\) of non-supports of Vj on C(xi,Xi+i) 
which satisfy B have to be scanned, in order to determine for 
each set a value w G Sj minimizing s(xj+i, w) and a value 
w' G -<Sj minimizing s(xi+i,w') + 1. This process takes 
0(\D(xi+i)\) for each value, leading to 0(d 2 ) for the whole 
domain. Since all the variables need to be scanned and for all 
the values in domains the quantities are stored, phase © takes 
0(nd 2 ) in time, and O(£d0 in space. 

Phases © and © take 0(EdO time each since all the 
domains have to be scanned. By Proposition 4, all the 
non-GAC values have been removed after this last phase. 

If B G {<,>,<, >,true}, -< is known. The worst-case 
time and space results come from Phase ©. The bottleneck 
stems from the fact that, when a domain D(x{ ) is scanned, the 
minimum and maximum number of C-stretches of each value 
are computed from scratch, while an incremental computa- 
tion would avoid to scan D(xi+i) for each value in D(xi). 
This observation leads to Property 6. Again, we focus on the 
minimum number of C-stretches on suffixes. Other cases are 
symmetrical. 

Notation 2. Given seq_bin(JV, X, C, B), x% G X, < i < n 
and a value Vj G D{xi), if 'i < n — 1, let Vj denote the set of 



integer values such that a value s(vj ,w) G Vj corresponds to 
each w G D(xi+\) and is equal to: 

• s(xi + i,w) if and only if w G Sj 

• s(xi+i,w) + 1 if and only if w G ~^Sj 

Within notation 2, the set Vj corresponds to the minimum 
number of stretches of values in D(xi+i) increased by one if 
they are non supports of value Vj with respect to C. 

Property 6. Given a counting-continuous instance 
of seq.bin(N, X, C, B) such that B G {<,>,<,>, true} 
and Xi € X, < i < n — 1, if the computation of 
rcan w ^u(x i+ i){ s { v ji w )) f or a ^ v j G D(xi) can be per- 
formed in 0(\D(xi+i)\) time then GAC can be achieved on 
seq_bin in O(SDi) time and space complexity. 

Proof. Applying Lemma 1 to the whole domain D(xi) 
takes 0(\D(xi+i)\) time. Storing the minimum number of 
stretches for each value in D(xi) requires 0(\D(xi)\) space. 
Phase © takes 0(£dO space and O(Sdi) time. □ 

When they are represented by a counting-continuous in- 
stance of seq_bin, the practical constraints mentioned in the 
introduction satisfy a condition that entails Property 6: Given 
Xi, it is possible to compute in 0(|£>(a;,-_|_i)|) the quantity 
m ^ n weD{x i+1 ) i s ( v Oi w )) f° r a first value i> G D(xi) and 
then, following the natural order of integers, to derive with 
a constant or amortized time complexity the quantity for the 
next value v\, and then the quantity for the next value vi, 
and so on. Thus, to obtain GAC in O(Edi) for all these 
constraints, we specialize Phase © in order to exploit such 
a property. We now detail how to proceed. 

With respect to the constraints mentioned in the intro- 
duction corresponding to instances of seq_bin which are 
not counting-continuous, the same time complexity can be 
reached but the algorithm does not enforces GAC. 

Thus, when seq_bin represents Change, Smooth or In- 
CREASING_Nvalue, computing mm we rj^ x . +1 ^(s(vo,w)) for 
the minimum value vq = mm(D(xi)) (respectively the max- 
imum value) can be performed by scanning the minimum 
number of C-stretches of values in 

We now study for Change, Smooth and Increas- 
inG-Nvalue how to efficiently compute the value 
™n weD{x . +l) (s(v k ,w)) of v k G D(xi), either directly 
or from the previous value imri weD r x . +1 ^s(vk-i, w)), in 
order to compute mm we r)( x ( s {vj , w)) for all Vj G D(xi) 
in 0(\D(xi)\) time and therefore achieve Phase © in 
0(S Di ). 

The Change constraint 

Section 3 showed a reformulation of Change(A^, X, CTR) 
as seq_bin(AT', X, C,true) A [N = N' - 1], where C is the 
opposite of CTR. 

— If C is ' > ' (the principle is similar for '<','>' and ' <'), 
the instance of seq_bin is counting-continuous, because B is 
true and C is monotonic. The monotonicity of C, with its 
corresponding total ordering -<, guarantees that given three 
consecutive variables Xi-i, Xi, Xj+i and v\ G D(xi-\),V2 G 
D(xi), «3 G D(xi+i), if (vi, V2) and (1*2, V3) both violate C, 
then we necessarily have v\ >- V2 >- V3. Therefore, changing 



value V2 by a new value v 2 such that v\ -< v' 2 (to satisfy C) 
entails v 2 >- v 2 , and thus still v 2 y V3 (which violates C). 
It is not possible to remove (or, symmetrically, to add) two 
violations of C within an instantiation only by changing the 
value of one variable. The instance of seq_bin is counting- 
continuous and thus the algorithm enforces GAC (by Propo- 
sition 4). 

To achieve step 3. in 0(D(xi)), we introduce 
two quantities lt(vj,Xi+i) and geq(vj,Xi+i) respec- 
tively equal to ram w£[min{D{Xt))tVj[ (s(x l+1 ,w)) and 
mm we [ v . imax ( D ( x .^(s(x i+ i, to)). The computation is 
performed in three steps: 

1. Starting from vq = min(.D(xj)), that is, the value hav- 
ing the smallest number of supports for C on Xi+i, 
compute lt(vj,Xi + i) in increasing order of Vj. Tak- 
ing advantage that, given a value Vj-i G D{xi) and 
the next value Vj G D(xi), [min(.D(aij)), is in- 
cluded in [min(D(xi)), Vj[. Therefore, the computation 
of all mm w& [ min ( D ( Xt )y Vj [(s{x t+ i,w)) can be amor- 
tized over D(xi+i). The time complexity for comput- 
ing U(vj, Xj+i) for all Vj G D(xi) is in 0(\D{xi)\ + 
\D(x i+1 )\). 

2. Similarly starting from vo = m&x(D(xi)), compute in- 
crementally geq(vj, ie-s+i) in decreasing order of Vj, in 

0(\D( Xi )\ + \D(x i+1 )\). 

3. Finally, for each Vj G D{xi), min w <£D(x i+1 )( s ( v j , w)) 
is equal to min(Zf(Vj-, 2^+1), geq(vj, Xi+i) + 1). 

Since step 3. takes 0(D(xi)), we get an overall time com- 
plexity for Phase © in 0(E D i). 

— If C is '=', '7^', or \xi — Xi + i\ < est (the latter case 
corresponds to the Smooth constraint), the filtering algo- 
rithm does not guarantees GAC because the corresponding 
instances of seq_bin are not counting-continuous. Step 3. can 
also be performed in 0(D(xi)), leading to an overall time 
complexity for Phase @ in 0(Y>Di): 

• If C is '=' then for each Vj G D(xi) there 
is a unique potential support for C on Xi+i, the 
value Vj. Therefore, by memorizing once the value 
vmini in D(xi+i) which corresponds to the small- 
est minimum numbers of C-stretches on the suf- 
fix starting at x i+1 : VV,-, min we£ >( a . +1 )(s(u,-,ui)) 
= min(s(.Ti + i, Vj), s(x.j+i, vmini) + 1), assuming 
s(xi+i, Vj) = +00 when Vj ^ D(xi+i). 

• If C is '7^' then for each Vj G D(xi) there is 
a single non support. By memorizing the two val- 
ues vmini and vmin 2 which minimize the mini- 
mum numbers of C-stretches on the suffix starting at 

Xj.fi, for any value Vj mm we D(x i+1 ) {s(vj ,w)) is equal 
to: min(s(ii+i,«mmi) + 1, s(xi+i, vmin2)) when 
vmini = Vj, and s(xi+i, vmini) otherwise. 

• Smooth is a variant of Change(./V, X, CTR), where 
Xi CTRxi+i is \xi — xi+i\ > est, est G N, that can 
be reformulated as SEQ_BiN(iV', X, C,true) A [N = 
N' — 1], where C is \xi — Xi+i\ < est. Assume 
Vq = mm(D(xi)) and we scan values in increasing or- 
der. Supports of values in D{xi) for |x, — Xj+i| < est 



define a set of sliding windows for which both the starts 
and the ends are increasing sequences (not necessarily 
strictly). Thus, mm^gs^. (s(vj, w)) can be computed for 
all Vj in D{xi) in 0(\D(xi)\) thanks to the ascending 
minima algorithm} Given a value Vj G D(xi) the set 
-iSj of non supports of Vj on \xi — a?;+i| < est is par- 
titioned in two sequences of values: a first sequence be- 
fore the smallest support and a second sequence after the 
largest support. While scanning values in D(xi) these 
two sequences correspond also to sliding windows on 
which the ascending minima algorithm can also be used. 

The Increasing JSf value constraint 

It is represented by seq_bin(A^, X, =, <), which is counting- 
continuous (see [Beldiceanu et al, 2010b] for more 
details). The algorithm enforces GAC. Since B is not 
true, we have to take into account B when evaluating 
m i n tuer>(x i+ i)( s (ji w )) f° r eacn v j & D(xi). Fortunately, 
we can start from v^ = max(D(xj)) and consider the 
decreasing order since B is '<'. In this case the set of 
supports on B can only increase as we scan D{xi). C is 
'=', then for each Vj G D(xi) there is a unique poten- 
tial support for C on Xi+\, the value Vj. We memorize 
once the value vmini in D(xi+i) which corresponds 
to the smallest minimum numbers of C-stretches on the 
suffix starting at x%+i, only on supports of the current 
value Vj G D(xi) on B. \/vj, iTim weD ^ x . +1 ^(s(vj,w)) 
= min(s(xi+i,Vj), s(xi+i, vmini) + 1), assuming 
s(xi + i,Vj) = +00 when Vj £ D(x i+ i). Since the set 
of supports on B only increases, vmini can be updated for 
each new value in D(xi) in O(l). 

6 Related Work 

Using automata, Change can be represented either by Reg- 
ular [Pesant, 2004] or by Cost-Regular [Demassey et 
al, 2006]. In the first case this leads to a GAC algo- 
rithm in 0(n 2 d 2 ) time [Beldiceanu et al., 2010a, pages 584- 
585, 1544-1545] (where d denotes the maximum domain 
size). In the second case the filtering algorithm of Cost- 
Regular does not achieve GAC. 

Bessiere et al. [Bessiere et al., 2008] presented an encod- 
ing of the CardPath constraint with Slide 2 . A similar re- 
formulation can be used for encoding SEQ_BiN(iV, X, C, B). 
Recall that SLiDEj(C, [xq, x-y, . . . , a;„_i]) holds if and only 
if C(xij, . . . , Xij+k-i) holds for < i < —f 1 - Following a 
schema similar to the one proposed in Section 4 of Bessiere et 
al. paper, SEQ_BiN(iV, X, C, B) can be represented by adding 
a variable N' and n variables [M , . . . , M n _i], with M = 
and AI n -i = N'. seq_bin(A^, X, C, B) is then reformlated 
by SLiDE 2 (C',[Af ,x ,Mi,xi,...,M n _i,x„_i]) A [N 1 = 
N - 1], where C = [-iC(xj, Xi + i) A B(xi, x i+ i) A M l+ i = 
Mi + 1] V [C(xi,x i+ i) A B( Xl ,x l+ i) A M l+ i = Mi]. Ac- 
cording to Section 6 of Bessiere et al. paper, GAC can be 
achieved thanks to a reformulation of SLIDE2, provided a 
complete propagation is performed on C, which is the case 
because B(xi, Xi + i) and C(x{, xi + i) involve the same vari- 
ables. The reformulation requires n additional intersection 

'See http://home.tiac.net/~cri/2001/slidingmin.html 



variables (one by sub-sequence [Mi,Xi]), on which 0(n) 
compatibility constraints between pairs of intersection vari- 
ables and 0(n) functional channelling constraints should 
hold. Arity of C is k = 4 and j = 2: the domain of 
an intersection variable contains 0(d k " : ') = 0(d 2 ) values 
(corresponding to binary tuples), where d is the maximum 
size of a domain. Enforcing GAC on a compatibility con- 
straint takes 0(d 3 ) time, while functional channelling con- 
straint take 0(d 2 ), leading to an overall time complexity 
0(nd 3 ) for enforcing arc-consistancy on the reformulation, 
corresponding to GAC for seq_bin. To compare such a time 
complexity 0(nd 3 ) with our algorithm, note that O(SDi) is 
upper-bounded by 0(nd). 

At last, some ad hoc techniques can be compared to our 
generic GAC algorithm, e.g., a GAC algorithm in 0{n 3 m) 
for Change [Hellsten, 2004, page 57], where m is the to- 
tal number of values in the domains of X. Moreover, 
the GAC algorithm for seq_bin generalizes to a class of 
counting constraints the ad-hoc GAC algorithm for Increas- 
ing_N value [Beldiceanu et ai, 2010b] without degrading 
time and space complexity in the case where seq_bin rep- 
resents Increasing JSf value. 

7 Conclusion 

Our contribution is a structural characterization of a class of 
counting constraints for which we come up with a general 
polytime filtering algorithm achieving GAC under some con- 
ditions, and a characterization of the property which makes 
such an algorithm linear in the sum of domain sizes. A still 
open question is whether it would be possible or not to ex- 
tend this class (e.g., considering n-ary constraints for B and 
C) without degrading complexity or giving up on GAC, in 
order to capture more constraints. 
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